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INVESTIGATION OF THE CONTAMINATION OF THE GOULD (2003) HALO SAMPLE 

Andrew Gould' 



o 
o 
o 

< 

o 



> 

cn 
od 

o 
o 



ABSTRACT 

A recent astroph posting argued that the lGouldl (|2003a|) halo sample is substantially contaminated with thick- 
disk stars, which would then "wash out" any signature of gr anularity in the halo velocity distribution due to 
streams. If correct, this would imply that the limits placed by iGouldl (l2003bl) on streams are not valid. Here I 
investigate such contamination using six different indicators, 

1) morphology of the underlying reduced proper motion diagram used to select halo stars, 

2) comparison of kinematic and parallax-based distance scales 

3) comparison of derived halo parameters for the Gould (2003al) sample with other determinations 

4) a precision color-color diagram for a random subsample 

5) the 3-dimensional velocity distribution of a random subsample 

6) metallicity distribution versus kinematic cuts on a random subsample 

I estimate that the contamination is of order 2 percent. Thus, the upper limits on the density of nearby streams 
derived by Gould (2003b) remain valid. In particular, at 95% confidence, no more than 5% of local halo stars 
(within about 300pc) are in any one coherent stream. Determining whether or not this local measurement is 
consistent with CDM remains an outstanding question. 

Subject headings: stars:halo - galaxies:substructure 



1. INTRODUCTION 

A critical test of the current picture of hierarchical struc- 
ture formation, is measuring the "granularity" or substruc- 
ture w it hin the halo as a function of Galactocentric distance. 
iGouldl (l2003bh determined that, contrary to naive expecta- 
tion within this framework, the local Galactic halo is remark- 
ably smooth, with only ^5% possible contribution from sub- 
structures. However, this conclusion is only as good as the 
underlving lGouldl (l2003ai) sample: significant contamination 
by thick-di sk stars would tend to "wash out" any substruc- 
ture signal (iKepley et al.ll2007h . Prompted by this concern, I 
demonstrate below through six different tests, that this sam- 
ple is not in fact seriously contaminated and therefore that the 
current (still valid) limit of ^ 5% substructure must be under- 
stood by current models of Galaxy formation. 

2. MORPHOLOGY OF THE RNLTT REDUCED PROPER 
MOTION DIAGRAM 

The lGouldl (l2003al) halo sample is derived from the revised 
New L uyte n Two-Tenths (rNLTT) catalog of iGould & SalimI 
0)0^ and 'Salim & GouM (l2003h . The argument made by 
Keplev et al. (2007) for contamination of this sample is that 
Ryan & N orrisI (11991 ') found such contamination among the 
halo ca r ididates th at they had extracted from the underlying 
iLuvtenI (11979. 19 80) NLTT catalog. 

Figures □ and |2] (both taken from I SaUm & Gouldl l2002h 
show the reduced proper motion (RPM) diagrams for the 
NLTT and rNLTT, respectively. Halo stars are clearly sep- 
arated from disk/thick-disk sta rs in rN LTT (used by Goulg 
|2003a), but not NLTT (used by lRyan & N orris 1991) . 

Figures [3] and |4] further confirm the difficulty of extract- 
ing a clean sample of halo stars from NLTT. They show the 
distributions of main-sequence stars and subdwarfs on the 
NLTT RPM. Clearly, there is no way to select subdwarfs from 
Luyten's original NLTT without substantial main-sequence 
contamination. 

' Department of Astronomy, Ohio State University, 140 W. 18th Ave., 
Columbus, OH 43210, USA; gould@astronomy.ohio-state.edu 
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(B-R)l 

Fig. 1 . — Reduced proper motion diagram for the original NLTT catalog, 
i.e., using the original proper motions and, more importantly, the original 
photographic (B, R) photometry. (Original B — R was sometimes given to 
only 1 decimal place: for these, s mall random numbers have been added to 
the color to permit display). From lSalim & Gouldl j2002l) . 



Subsequently, ISalim & Gouldl (l2003h intro duced a further 
refinement of the optical/infrared RPM of Sali m & Gouldl 
(I2OO2I) by defining a discriminator rj that depends on both 
RPM and Galactic latitude b: 



77 = y + 51og^-3.1(y-y)-1.47|sin/7|-2.73 



(1) 



where fi is the proper motion in arcsec per year. They identi- 
fied halo stars as dominant in the range 

0<?7<5.15 [Salim & Gould (2003) : halo stars] (2) 

Figure |5] adapted from ISalim & Gouldl ( 120031) shows (in 
black) the distribution of r] in the color range 2.25 <V -J < 



15 



V-J 

Fig. 2. — Reduced proper motion diagr am for revised rNLTT catalog 
IGould & SaIin]|2003t[S"alim & Gouldll2003!) i.e., u.sing the new proper mo- 
tions and, more importantly, new 2MASS J photometry to enable a muc h 
broader-baseline and more accurate V-J color From iSalim & Gould )200A . 
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Fig. 3. — Reduced proper motion diagram for the original NLTT catalog. 
I.e., same as Fig.Hrestricted to rNLTTHNLTT, but with main-sequence stars 
(as determined from Fig.|2) shown in red. From iSalim & Gould Ii2002i^ . 



3.25. This confirms quantitatively the visual impression from 
Figure |2] that there is a clear valley between the subdwarfs 
and the main-sequence stars. To further clarify the situation 
I have estimated the subdwarf/main-sequence breakdown in 
the region of overlap as follows (green and red histograms in 
figure). First, at the valley minimum (77 = 0) I assigned half 
the stars to each population. Second, for 77 < 0, I estimated 
the falling profile of subdwarfs to be the mirror image of their 
rising profile at high 77. Third, for 77 > 0, 1 estimated the falling 
profile of main-sequence stars to be the mirror image their ris- 
ing profile at low 7/. 

In or der to obtain a pure sample of halo stars, iGouldl 
( l2003ah further restricted the range of rj relative to equation 



(B-R)l 

Fig. 4. — Reduced proper motion diagram for the original NLTT catalog, 
i.e., same as Fig.rH'estricted to rNLTTHNLTT, but with subdwarfs (as deter- 
mined from Fig.|2) shown in green. By comparing this figure with Fig. [3] it 
is clear that it is impossible to design selection criteria that would recover a 
large number of subdwarfs from the original NLTT without substantial con- 
tamination by main-sequence stars. From .Salim & Gould l .2002i) . 



dU, using 

1 <77<4.15 



[Gould (2003) : secure halo stars] (3) 



to avoid main-sequence stars at high rj and white dwarfs at low 
7]. From Figure|5] one may estimate that this choice generates 
roughly 2% contamination. 

While I will give several other independent arguments that 
the sample is not seriously contaminated, this one is the 
strongest and most quantitative. 

3. COMPARISON OF PARALLAX AND KINEMATIC 
DISTANCE SCALES 

A second check comes by comparing th e kinematic-based 
distance scale derived by iGouldl (l2003al) with trigonomet- 
ric par allaxes found in the literature. To establish dis- 
tances, IGouldl ( l2003al) fit for the two parameters of a lin- 
ear color-magnitude relation, while enforcing a mean mo- 
tion of halo stars relative to Sun of U2 = -2I6.6kms~', so as 
to force agreement with the determination of this parameter 
by G ould & Popowski (1998) based on halo RR Lyrae stars. 
IGouldl (l2003al) then compared this with the color-magnitude 
relati on derived from hal o star s with trigon ometric parallaxes 
from lMonet etal.1 d 1 9921) and iGizisI (Il997h . As can be seen 
from Figure |6] these relations are virtually identical. If the 
sample were contaminated with thick-disk stars (whose mean 
motion U2 motion relative to the Sun is of order 5 times 
smaller than halo stars) by even 10%, then this would cause 
an error in the distance scale of 0.80 x 10% = 8%, yielding an 
offset between the parallax and kinematic relations of about 
0.17 mag. The offset is clearly much smaller than this. 

4. COMPARISON OF HALO KINEMATICS 

It is not only the asymmetric drift o f the halo that is cor- 
rectly reproduced by the Gould (2003a) analysis, but also the 
velocity dispersions. Gould (2003a ) found dispersions in the 
radial, rotation, and vertical directions of 168 ± 1, 113 ±2, and 
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Fig. 5. — Distribution of "population discriminator" r] for rNLTT stars 
with 2.25 < V —J < 3.25. Tliere are two clear peaks corresponding to main- 
sequence stars (red) and subdwaii's (green). In the overlap the bins have been 
divided by symmetrizing (see text). Adapted from Salim & Gouldi >2003i) . 
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Fig 6.— CMD of stars with trigonometric parallaxes ('Monet et al."192 
IGizis|[T997l) . separated into various classes using the Salim & Gould (200; 
RPM-discriminator rj. The solid line shows the best fit to the halo stars in 
the data while the dashed line shows the fit based on kinematic analysis of 
thelGould (2003a) sample. They are virtually identical. If there were serious 
contamination by thick disk stars, the dashed line would have been driven to 
brighter mags. From , Gould C2003a) 
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Fig. 7. — BVI color-color plot of 564 stars selected from rNLTT using the 
same 1 < < 4. 15 criterion apphed by Gould (2003a). The sample is metal 
poor and relatively homogeneous. From , Marshall (.2007i) . 

5. PRECISION COLOR-COLOR DIAGRAM OF A 
RANDOM SUBSAMPLE 

Marshall has obtained precise photometric data for 564 
stars selected by applying the same 1 < 77 < 4.15 criterion 
used by Gould (2003a). Her Figure 9 (reproduced here as 
Fig-ID shows a very tight color-color relation for these stars, 
consistent with a metal-poor population. It is clear that the 
stars extracted by this criterion are not very heterogeneous, as 
would have been expected were they seriously contaminated. 

6. 3-D VELOCITY DISTRIBUTION OF A RANDOM 
SUBSAMPLE 

J. Marshall (private communication 2007) has also obtained 
radial velocities for 295 of the stars in her sample, which (to- 
gether with her excellent photometry and the rNLTT proper 
motions) permit her to make 3 -dimensional velocity estimates 
for each star. The results are shown in Figure [8] The distri- 
bution is not significantly contaminated by thick-disk stars, 
which would appear as an overdensity centered at (U,V,W) ^ 
(0,190,0)kms-i. 



89 ± 2 km s , respectively. These values are quite compatible 
with other determinations. For example Gould & PopowsQ 
(0.998) found 171 ± 10, 99 ± 8, and 90 ± 7 from a much 
smaller sample of halo RR Lyrae stars. Serious contamina- 
tion by the kinematically much cooler population of thick- 
disk stars would have tended to drive down these dispersions 
(assuming the distance scale remained in agreement with the 
parallax stars of Fig.|6]l. 
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Fig. 8. — UVW velocities (in Galactocentric frame) of 295 stai's selected 
from rNLTT using the same 1 < 77 < 4. 15 criterion adopted by Gould ( 2003;j). 
There is no significant contamination by thick-disk stars. Kindly provided by 
J. Marshall in advance of publication. 
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Fig. 9. — Metallicity distribution for 239 stars selected from rNLTT us- 
ing the same 1 < r; < 4. 15 criterion adopted by Gould (2003a) (black) and 
for a "kinematically secure halo" subsample of 194 stars (red). The latter 
histogram is scaled up by 239/194 so that the sum of the bins is the same. 
The two histograms are essentially identical, while the analogous Figure 2 
of Ryan & Norris 1 1991) shows that their NLTT-derived sample was heavily 
contaminated by thick-disk stars. Kindly provided by J. Marshall in advance 
of publication. 



which they showed that the halo candidates that they had 
selected from NLTT were heavily contaminated by thick- 
disk stars. Figure [9l is the direct analog of Figure 2 from 
iRvan & NorrisI ( Il991h . It shows the metalHcity for all 239 
stars and also for a subsample of 194 "kinematically se- 
cure halo stars" defined as the union of stars with either 
vj^ > 220 km s"' or velocity in the rotation direction V < 
-220kms~'. The first point to note is that 82% of the stars 
are "secure halo" by this definition. The second point is that 
the metallicity distributio ns of the two s amples are essentially 
identical. By contrast, Ryan & NorrisI ('1991) found a huge 
tail of higher-metallicity stars in their full sample, which dis- 
appeared when they implemented the same "secure halo star" 
cut. This brings the argument full circle. 



8. CONCLUSION 

Contamination of the original iGouldl ( l2003ah halo sample 
by thick-disk stars is likely to be of order 2%. This implies 
that the effect of streams is not "washed out" by thick-disk 
contamination . Hence, the limits on granularity found by 
iGouldl (l2003bl) in this sample imply corresponding limits on 
streams in the Galactic halo. In particular, at 95% confidence, 
no more than 5% of local halo stars (within about 300pc) are 
in any one coherent stream. 



7. METALLICITY DISTRIBUTION VS. KINEMATIC 
CUTS 

J. Marshall (private communication 2007) was able to make 
preliminary metallicity measurements for 239 of the 295 halo 
candidates mentioned in § |6] Th ese permit a direct com par- 
ison with the test performed by iRvan & NorrisI (1199 ll) . by 
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